Giant Microwave Absorption in Metallic Grains: Relaxation Mechanism 
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We show that the low frequency microwave absorption of an ensemble of small metallic grains at 
low temperatures is dominated by a mesoscopic relaxation mechanism. Giant positive magnetore- 
sistance and very strong temperature dependence of the microwave conductivity is predicted. 
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Microwave absorption in an ensemble of metallic grains 
has been investigated both theoretically and experimen- 
tally in many papers. (See for Refs. At high 
temperatures quantum effects can be neglected and the 
effective microwave conductivity at frequency lo is given 
by the classical Debye formula, 
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lo 



(1) 



Here a c i is the classical Drude conductivity determined 
by elastic scattering. At low temperatures the quantum 
nature of electronic states in grains becomes essential and 
statistics of the electron levels determines the microwave 
absorption [|l]-|4|. However all quantum effects consid- 
ered so far result in the microwave conductivity of the 
order or smaller than <tjj. In this Letter we discuss the 
mechanism of microwave absorption which can be much 
stronger than the classical one given by Eq. (1). 

At sufficiently high radiation frequency lo the absorp- 
tion is dominated by the resonant mechanism, i.e., by 
the direct transitions between electron levels and 
the microwave conductivity is determined by the proba- 
bility density P(s) of energy spacing s between adjacent 
levels. This quantity is usually described by the Wigner- 
Dyson distribution. At s < A, where A is the mean 
level spacing, the probability density behaves as § (see 
the insert on Fig. 1) 



P(s) = Cps^/A^ 



(2) 



The exponent [3 is determined by the global symme- 
try of the system. For the orthogonal ensemble (3=1 
and Cp = 7r 2 /6. In the unitary case, (3 = 2 and 
C/3 = 7r 2 /3. Finally, the symplectic ensemble is char- 
acterized by (3 = 4 and Cp = 167r 2 /135. 

It is well known, however, that absorption of radiation 
in two level systems is determined at low lo by the relax- 
ation mechanism. Mandelstam and Leontovich proposed 
this mechanism for sound attenuation in gases with slow 
internal degrees of freedom For microwave absorp- 
tion in doped semiconductors relaxation absorption was 
suggested by Pollack and Geballe [Q . 

The contribution of this mechanism to the microwave 
conductivity of metallic grains or depends strongly on 



the relation between A and the level broadening T. For 
r > A, the relaxation mechanism of absorption in gran- 
ular metals was shown in Ref. Q to determine under 
certain conditions both the microwave conductivity and 
the magnetoconductance of metallic grains. Here we con- 
sider the opposite limit r <C A when the energy levels 
are well resolved. It turns out that the microwave ab- 
sorption ctr as well as the magnetoconductance are giant 
in this limit, i.e., an exceeds both Eq. (1) and results of 
Ref. d) by several orders of magnitude. 

We now discuss the physical picture of the relaxation 
absorption of applied electric field of frequency lo. In 
adiabatic approximation temporal energy levels tiit) os- 
cillates with the same frequency. Since the populations 
of the energy levels follows adiabatically the motion of 
the levels themselves, the electron distribution becomes 
nonequilibrium without any interlevel electron transi- 
tions. Relaxation of the nonequilibrium distribution due 
to inelastic processes leads to the entropy production and 
therefore to the absorption of the energy of the external 
field. 

To simplify the situation, we will restrict ourselves to 
the low temperature case when T <^ A. (Later we will 
also discuss the situation at T <C A < T qualitatively.) 
We will see that at lo <C T the relaxation mechanism of 
absorption is determined by rare grains where the first 
excited level is separated from the ground state by small 
energy s <C T. This implies that the ground and the first 
excited states in the grains form two level systems which 
are effectively isolated from the rest of the spectrum. 

The population N(t) of the first excited state with the 
energy s(t) is governed by the equation 



dN(t) = N(t) - N (s(t)) 
dt T f 



(3) 



Here No(s(t)) = (l+exp(s(i)/T)) _1 is the adiabatic equi- 
librium population of the first excited level and r e is the 
relaxation time of the level due to the electron-phonon 
interaction. Eq. (^) is applicable provided s>T. 

It is crucial for our discussion that the relaxation rate 
for a two level system vanishes with s. For a two level 
system embedded into 3D insulating environment with 
phonons, r e ~ s -4 for s — * 0. The fact that t £ (s) diverges 
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so rapidly for small s leads, as we will see, to divergences 
in <jr. To provide the proper cutoff we need to include 
the contribution to r" 1 from the processes that involve 
the second excited state. The latter is typically sepa- 
rated in energy by A > T and we get a finite, though 
exponentially small, relaxation rate at s = 0: 



T £ (s) T A T m 



where 
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flD is the Debye frequency in the metal, v s denotes 
the velocity of sound, L is the size of the sample 
and I is the elastic mean free path. To interpret 
Eq. (||) we present the first term in Eq. ([|) in the form 
(s 3 /n 2 D )(A/s)(kl){kL) 2 (T/s) 1 where k = {s/hv s ) is the 
wave number of a phonon with the energy s. The first 
factor here corresponds to the conventional expression 
for the relaxation rate in clean bulk metal at ~ T. 
We are considering a transition between two particular 
levels without summing over the final states (as is usual 
in the bulk case); this gives rise to the second factor. In 
the dirty case {kl < 1) the electron phonon relaxation 
time is known to be suppressed H (third factor). If the 
phonon wave length exceeds the system size L as well, the 
electron-phonon coupling in the dipole approximation is 
reduced, and the fourth factor appears. Finally at s < T 
the large population number of phonons leads to the last 
factor. 

To evaluate <tr we consider the power Q absorbed in 
a grain. Given the amplitude of the microwave electric 
field E and the volume of the grain V, 
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E 2 V 



(0) 



In the ohmic approximation we obtain the Debye-type 
expression for the averaged microwave conductivity, 

dsP(s)^-^^. (7) 



or 
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ds 1 + (cjt £ (s)) 2 



Here (...) stands for the averaging over the random scat- 
tering potential. The random matrix elements were as- 
sumed to be statistically independent from the spectrum. 
This assumption can be verified by a straightforward cal- 
culation. We will present this calculation elsewhere. 

The distribution P(s) of the nearest neighbor spacing 
s is known to be determined by the global (spin and/or 
time reversal) symmetry of the system. a R , because of 
its dependence on P(s) as seen in Eq. ([?]), should be very 
sensitive to weak magnetic field H and to spin-orbit scat- 
tering rate 1/t so . Here we restrict ourselves only to the 
three asymptotic regimes: (i) orthogonal, when H = 



and t so — oo, (ii) unitary, when the magnetic field is 
strong and (hi) symplectic (H — 0, while t so is short.) 
We also discuss the particularly interesting case of the 
combination of strong spin-orbit scattering and a weak 
magnetic field. With this exception we do not present 
here explicit formulae for the absorption in the crossover 
regimes. Here we will consider only linear absorption. 
Substitution of Eq. (||) into Eq. (^) gives 
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where 
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For orthogonal (/3=l)and unitary (/3=2) ensembles of 
metallic grains (t so s* 1) the typical value of s deter- 
mining the integral contribution in Eq. (7) is 



s - s* = A(max[wTo,To/T m ]) 1/4 , 
provided s* <c T, A. In this limit, or is 



(10) 



= ^0 // d£.)2v/Zm jl/2/3 / ( arctan ( WTm )) (11) 

WVT^dE' A t ; J/n v ;; v ; 

with 

Mk) ^cosec^HcosAO^cos^). (12) 

The amplitude of electric field decays over Thomas- 
Fermi screening length vq = D finer c i L from its 
vacuum value E down to Euj/a c i <C E. (Here D is the 
diffusion constant of electrons.) This small electric field 
in the bulk gives the main contribution to the classical 
Debye formula Eq. (1). On the contrary, the sensitivity 
Eq. (12) is determined by the small region of the width of 
order of tq near the surface, where the value of the elec- 
tric field is of the order of E. A calculation analogous to 
that in Ref. B gives 



ds 



4^4 



e r. 



y dE' ' (3a cl V 
As a result, in the orthogonal case 



o n I v a > i- 

a R = Cacr D x< l /_i_u/ 



(J) 1 / 2 exp(A) W T m «l 

M — ) V2 U>T m » 1 



(13) 



. (14) 



while for the unitary ensemble, 
a R Caa D x j 2 )( _^ )1/4 ^ >> 1 , (15) 



where 



a = 7rA 2 T /372T 



(16) 
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and C is a constant of order unity which depends on the 
geometry of the grain and on the direction of the electric 
field. If the grain is cubic and the microwave field is per- 
pendicular to its face, then C = 1. Therefore in both the 
orthogonal and unitary cases ctr is much larger than the 
classical Debye conductivity Eq. (1). 

In the absence of electron-electron Coulomb interac- 
tion, a similar approach is applicable even at T > A > T 
when many well resolved levels participate the absorp- 
tion. In this case <jr ~ a'ajj and a' = e TqV T e (T). 
This matches the results || at A ~ r. 

As usual, one can drive an orthogonal system into a 
unitary one by applying magnetic field H. This leads to 
a giant negative magnetoconductance. Here we consider 
it only qualitatively and only at wtm -C 1. The field 
needed to reduce from a° R to a R is 



H ~ H(s*) 



(17) 



where given the cross-section of the grain A and its di- 
mensionless (in units e 2 /7i) conductance g, the charac- 
teristic field Ho is 



Hn = 



he 

eAgi 



(18) 



The field H(s*) reduces o\r dramatically from a° R to a R . 

Now let the spin-orbit scattering rate t." 1 exceed A. In 
this symplectic case P(s) ~ s 4 /A 5 and the integral over 
s in Eq. (6) is determined by s ~ T. At low u> <C t~ 
(which means s* <C T) we have 

a% = C^fa D a^; if w « (T/A)V- (19) 

At small frequencies a s R turns out to be temperature in- 
dependent and larger than <td by the factor aT/A ~ 
roA > 1. At high frequencies the condition s* <C T is 
violated. In the interval T < s* < A (or (T/A) 4 < ujt < 
1) a r is u>- independent for all ensembles and 



—^—(T/Ar^T 
16np(p + 5) r 



(20) 



In a symplectic grain each energy level is double degen- 
erated due to the T-invariance. An applied magnetic field 
splits this Kramers doublet and therefore increases P(s) 
at small s. According to Kravtsov and Zirnbauer jl(J 



P(s,H)-P(s,0) 



8(H)- 



■ exp 



s{Hf 



Therefore the distribution function has a peak at 

H 



s = s(H) = A 



Ho 



(21) 



(22) 



and P(s) ~ s 2 below this peak. Calculation of the sen- 
sitivity ((ds/dE) 2 ) and T e for a split Kramers doublet 



require some caution: without magnetic field both the 
external electric field and lattice deformation are unable 
to split a doublet (for example ds/dE = for H = 0). 
Therefore at H < H 



((J|) 2 >-e 2 r 2 -^(|-) 2 
dE Ana e i H 

1 



r £ tq A l iV 



(23) 
(24) 



Substituting Eq. ( ^lj ) into Eq. (0) we get the expres- 
sion for grain magnetoconductance in the limit wt ( <C 1 
in the presence of substantial spin-orbit scattering. For 
s(H) > s*, i.e., H 4 /H$ > wto,t /t to , we get 

a R (H) - a R (0) ~ (^) 3 ^ » <7 s (il = 0). (25) 

Thus we have a giant positive magnetoconductance that 
increases rapidly with the decreasing magnetic field until 
H > H s*/A. Therefore at H - H s*/A, the magne- 
toresistance has a sharp maximum. The magnetic field 
dependence of aft both for t so — -» oo and for short t so is 
qualitatively illustrated in Fig. 1. 

Let us now estimate (Tr. Consider a metallic grain 
with A ~ IK (the size is about L - 150A). To can be 
estimated as To ~ 10 -7 sec. Then we get a ~ 10 2 at 
T ~ 0.3K. Therefore, at lu ~ 10 5 Hz we predict cr^ in 
the orthogonal case to be about three orders of magni- 
tude larger than an- Under these conditions s* ~ 0.3K 
and # ~ 10T (at g ~ 10), i.e., H ~ 3T will make the 
absorption several times smaller. 

Note that the magnetoconductance is negative and 
positive in the orthogonal and symplectic cases respec- 
tively, i.e., the sign of the magnetoconductance is op- 
posite to usual weak localization |ll|,E3]. This is natural 
since classical microwave conductivity Eq. (1) is inversely 
proportional to u c \. 

So far we have neglected effects of the Coulomb interac- 
tion between electrons. First of all, our consideration was 
based on the Wigner-Dyson distribution P(s). Its appli- 
cability in the presence of interactions is not as clear as 
in their absence. Experimental study of the microwave 
absorption could provide information on this probability 
distribution. 

Another Coulomb effect to consider is the contribution 
of the electron-electron interaction to the level broaden- 
ing r. We argue that in the limit T <C A this contri- 
bution does not exist and T is determined entirely by 
the phonon scattering T — t~ . Coulomb interactions 
can only renormalize the one-electron energy of the first 
excited state but not to broaden it, because the inelastic 
transitions due to electron-electron interaction are forbid- 
den by the energy conservation law. Therefore Eq. (|t]) 
remains valid even for strongly correlated electrons. 

As it has been mentioned above, without Coulomb in- 
teraction there is a big interval r ~ t" 1 < A < T where 
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our two level approximation is not valid and many levels 
in each grain contribute to the absorption. These levels 
are still well resolved. On the other hand the electron- 
electron interaction is known to dominate T at high tem- 
peratures: r ee > A for T > A pi. This makes cal- 
culations for A < T < A(T A) 1 '^anything but triv- 
ial. Indeed, calculated in Ref. T ee is the broadening 
of one-particle states. The exact many-electron states 
can be broadened only by phonons. However, the spac- 
ing between these exact states A ex (e) is of the order 
of A only for several low energy excitations. With in- 
creasing the energy e the spacing A ex (e) decreases very 
quickly. We find it possible that the results of (s) are 
valid until T > T* where T* is determined by the rela- 
tion A ex (T*) = T p h(T*). If so, the matching should take 
place in the interval A < T < T* . We have to admit that 
much better understanding of the relation between exact 
states and quasiparticles and of electron-phonon interac- 
tion in the closed system is needed in order to develop a 
theory of microwave absorption at A < T < A(toA) 1 / 5 . 

In addition to the absorption of the electric field, one 
can consider the effect of an a.c. magnetic field. Here we 
only mention results of Ref. jl3|. The authors of these 
papers assumed that the broadening of the energy levels 
is energy independent and found that the energy dissipa- 
tion in the time dependent magnetic field can be much 
larger than the classical one at A 3> V. The sign of the 
magnetoresistance in this case turns out to be the same 
as we found here, i.e., magnetoresistance is positive in an 
orthogonal ensemble and negative in a symplcctic one. 
We believe, however, that the energy dependence of T is 
crucial for this problem as well, and we plan to present 
a quantitative theory elsewhere. 
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FIG. 1. Magnetic field dependence of the microwave con- 
ductivity at (a) strong spin-orbit scattering, and (b) t 30 — > oo. 
Note a sharp and high peak in the first case due to the 
splitting of the Kramers doublets. Parameters cr'^ u ' s , s* , A 
and Ho are determined in the text. Insert: P(s) at small 
s (s <C A): (1) orthogonal, (2) unitary, (3) symplectic (4) 
t so A > 1 and finite H. 
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